UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PACE  (Whan  OM«  Bnimrmd) 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3.  RECIPIENT'S  CATALOG  NUMBER 


■^TOMATIC  TARGET  DESIGNATION  / 


7.  authorc*; 

1 RichEu:d  D.j: 


I Hilton 


•.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Naval  Surface  Weapons  Center  (K61) 
Dedilgren  Lal>oratory 

qren.  Virginia 22448 


II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Chief  of  Naval  Material 
Department  of  the  Navy 


•.  CONTRACT  OR  GRANT  NUMBERf*; 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  • WORK  UNIT  NUMBERS 


62766N; ZF61312 : ZF61312001 


monitoring  agency  name  a AODRESSflf  <M//*ranl  /root  Controlling  Otlleo)  IS.  SECURITY  CLASS,  (ol  thio  ropor/J 

UNCLASSIFIED 


16.  DISTRIBUTION  STATEMENT  (ol  thlo  Koport) 


6 i F'  I 


Approved  for  public  release;  distribution  unlimited 


Television  Guidsuice 
Automatic  Target  Designation 
Terrain  Models 

Computerized  Simulation  


RACT  fConflmi*  on  rovoroo  old*  It  noeoooor^  and  Identity  by  block  mmbof) 

Television  guidance  systems  have  traditionally  relied  upon  a human 
operator  to  designate  a target  before  weapon  release.  This  study  uses 
computer  simulations,  which  include  a simulated  terrain  model,  to  evolve 
on-boEu:d  hardware  techniques  for  automatically  designating  targets  after 
weapon  release.  Projectiles  that  use  128-line  tele'vlsion  guidemce  that 
automatically  designates  moving  objects  in  a Izmd  tzurget  area  are  successfully 


EDITION  OF  I NOV  61  IS  OBSOLETE 
S/N  0102-  LF- 014- 6601 


UNCLASSIFIED 

security  CLASSIFICATION  OF  THIS  PAGE 


('BW«i  OMa  Elilwam 


if>^U 


FOREWORD 


Automatic  Target  Designation  was  a FY  1978  lED  project  designed  to  illus- 
trate the  applicability  of  television  guidance  technology  to  a "fire-and- 
forget"  guided  projectile. 

This  report  was  reviewed  by  C.  T.  Shelton,  Head,  Computer  Engineering 
Branch  and  C.  J.  LaMonica,  Head,  Computer  Facilities  Division. 


Released  by: 

R.  A.  NIEMANN,  Head 
Strategic  Systems  Department 


iii 


CONTENTS 


AUTOMATIC  TARGET  DESIGNATION/MOVEMENT  DETECTION 
VTRAJ3 — A PROGRAM  TO  GENERATE  A VIDEO  SIGNAL 

TERRAIN  MODEL  

OPTICS  SIMULATION  

trajectory  SIMULATION  

MOVING  TARGET  

TARGET  DISCRIMINATION  

CONCLUSIONS  

APPENDIX  A — PROGRAM  LISTINGS  

DISTRIBUTION 


1 

2 

2 

9 

11 

13 


13 

29 


35 


fwcEDiNo  page  BUMC-NOT 


V 


ILLUSTRATIONS 


Figure 


Outline  of  the  Computer  Study  

Plan  Sketch  of  the  Terrain  Model  

Data  Card  Representing  a Portion  of  the  Terrain  Model  MAP  . 

Typical  Sketch  Used  to  Formulate  Subpictures  

Data  Card  Representing  the  Upper  Half  of  Subpicture  00.  . . 
Output  of  Program  DISPLAY  Showing  the  64  Subpictures  . . . 

Coordinates  Used  for  the  Optics  Simulation  

128^  Pixel  Data  Format  With  the  Corresponding  35mm  Film 

Format  

Trapezoidal  Projection  of  a Pixel  on  the  Terrain  Model  . . 

Typical  Trajectory  and  View  Angle  

Computer -Gene rated  Pictures  of  an  Area  Containing  a Moving 

Target  

Video  Processor  Simulation  

Sample  Magnification  Process  of  a 15^  Segment  Expanded  into 

a 16^  Segment  

Superimposed  Pixels  from  an  Input  Picture  and  a Magnified 

Output  Picture  

Hardware  Array  to  Effect  Image  Transformation  

Typical  Picture  (Top)  Operated  Upon  by  Various  Values  of 

Magnification  

Data  Flow  of  Simulation  

Flow  Chart  of  XMARK  

Correlation  Function  Curve  Pit  

Output  of  Program  XMARK  Showing  a Target  Designation  . . . 

Four  Designations  of  Moving  Targets  

Alternate  Hardware  Configuration  to  Eliminate  Image  Trans- 
formation Device  


SI  CONVERSION 


Multiply 


To  Obtain 


0.3048  meter 
1.6094  kilometer 
0.9144  meter 


HffiCEDlNO  PAGE  BUMC.N0T  FIl>iD 


AUTOMATIC  TARGET  DESIGNATION/MOVEMENT  DETECTION 


The  goal  of  the  automatic  target  designation  effort  is  to  overcome  the 
limitation  usually  imposed  by  television  guidance  systems  which  require  a 
human  operator  to  place  a curser  on  an  object  to  be  designated  as  a target. 

With  the  cursor  in  place,  television  guidance  systems  have  been  proven  suc- 
cessful in  guiding  to  a target  without  further  assistance.  A method  is  sought 
by  which  television  camera  output  may  be  processed  so  as  to  designate  a target 
without  human  intervention. 

The  process  we  seek,  however,  must  take  the  place  of  a human  brain  that 
is  making  judgments  based  upon  visual  perception.  The  quantity  of  data  avail- 
able from  a television  sensor  is  very  high  (typically  around  10®  bits  per 
second) , and  the  processing  is  quite  complex;  thus,  the  search  was  approached 
in  the  following  manner : 

1.  To  attenpt  to  process  data  so  as  to  discriminate  an  object  or  objects 
that  are  moving  with  respect  to  a complex  background  of  stationary  features 
(land  targets) 

2.  To  carry  out  the  search  in  the  form  of  computer  simulations  (avoid 

the  use  of  special  facilities  such  as  the  Army's  terrain  model  at  Huntsville, 
Alabama) 

3.  To  use  simplified  vehicle  motions  such  as  a simple  trajectory  with 
air  drag  but  without  vibrations  or  vehicle  roll 

4.  To  process  data  in  such  a way  as  to  simulate  configurations  of  hard- 
ware components  that  are  either  available  now  or  are  anticipated  in 

the  near  future 

5.  To  outline  broadly  any  design  features  that  may  not  be  standard 
practice 

Figure  1 illustrates  the  broad  outline  of  the  computer  simulation. 

The  blocks  labeled  trajectory,  optics,  target,  and  background  are  encompassed 
in  a single  program  named  VTRAJ3  which  has  the  objective  of  providing  simu- 
lated video  signals  from  which  target  designations  are  extracted.  With  a 
number  of  simulated  flights  on  disk  file,  other  programs  are  used  to  process 
and  edit  the  video  data.  The  data  reduce  and  printout  functions  are  accom- 
plished through  the  use  of  a simple  transcribing  program  and  the  Stromberg 
Carlson  4060  Cathode  Ray  Tube  Graphics  system.  All  output  is  rendered  on 
35mm  film  and  is  used  as  a projection  slide  or  is  printed  by  the  Photography 
Lab . 


A process  capable  of  designating  moving  targets  has  been  found.  The 
following  sections  describe  the  programs  and  processes  in  detail,  and  the 
final  section  evaluates  the  findings  in  quantitative  terms. 


1 


Figure  1.  Outline  of  the  Computer  Study 


VTRAJ3--A  PROGRAM  TO  GENERATE  A VIDEO  SIGNAL 

The  Fortran  program  designated  VTRAJ3  was  written  to  generate  files 
of  data  representing  the  video  signal  that  would  eminate  from  a 128-line 
TV  camera  mounted  in  a projectile  that  is  descending  along  a ballistic  path 
onto  a terrain  target  area.  The  camera  is  mounted  so  as  to  look  along  the 
velocity  vector;  and,  it  is  assumed  that  the  vehicle  is  roll  stabilized  with 
the  camera  upright.  There  is  a provision  for  superimposing  a moving  target 
upon  the  surface  of  the  terrain. 


TERRAIN  MODEL 

A terrain  model  was  required  to  present  the  moving  target  against  a 
realistic  background.  Since  no  physical  terrain  model  was  readily  available, 
a terrain  model  simulation  was  built  into  the  computer  program.  The  model 
consists  of  a digitized  mosaic  representing  a rectangular  area  1920  by  4800  yd. 
Each  square  yard  within  this  area  is  assigned  a singie  octal-digit  gray  level 
with  0 representing  black  and  1_  representing  white.  The  mosaic  of  gray  squares 
forms  an  isometric  drawing  of  a terrain  model  with  features  portrayed  as 
seen  from  above,  looking  down  at  a 45°  angle. 

In  use,  the  terrain  model  is  a table-look-up  subroutine  which  is  given 
integer  (I,  J)  coordinates  in  yards  and  returns  with  the  corresponding  gray 
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level.  To  avoid  excessive  memory  requirements,  the  table- look-up  is  completed 
in  two  steps. 

The  entire  terrain  model  area  is  represented  by  a two-dimensional  array 
called  MAP  (48,  120).  Each  element  in  MAP  represents  a square  area,  40  yd 
on  a side.  The  representation  is  a number  in  the  range  (0  to  77) g which 
designates  one  of  64  standard  subpictures.  Each  subpicture  is  a table  in 
which  each  of  the  1600  one-yard  squares  represented  is  a single  octal  digit. 
The  octal  digits  are  packed  20  to  a 60-bit  CDC  6700  word  so  that  the  entire 
subpicture  set  is  contained  in  an  array  called  NPIX  (80,  64) . 

Figure  2 is  a plan  sketch  of  the  terrain  model.  Each  of  the  graph  paper 
squares  shown  represents  a subpicture.  The  sketch  within  each  square  is 
a rough  rendering  of  one  of  the  64  standard  subpictures.  A corresponding 
set  of  data  cards  was  then  made  with  the  input  format  containing  24  array 
elements  per  card.  The  entire  array  was  punched  on  a deck  of  240  cards. 

Figure  3 shows  a section  of  the  data  preparation.  The  first  120  cards  con- 
tain data  for  the  left  half  of  the  terrain  model,  and  the  second  120  cards 
contain  data  for  the  right  half.  The  card  deck  is  given  the  permanent  file 
name  MAP  and  the  local  file  name  TAPE2.  The  file  is  read  under  control  of 
the  following  Fortran  statements: 

READ  (2,2) ( (IMAP(I,J) ,1=1,24) ,J=1,120) 

READ  (2,2) ( (IMAP(I,J) ,1=25,48) ,J=1,120) 

2 F0RMAT  (2402) 

The  subpictures  began  as  sketches  drawn  in  groups  of  four  on  graph  paper 
with  each  square  representing  1 yd^.  The  sketches  were  then  digitized  by 
sight  onto  keypunch  forms  and  were  made  into  a permanent  file  named  SUBPIX. 
Each  subpicture  occupies  40  cards.  Figure  4 is  a typical  sketch,  containing 
roads  and  houses,  of  four  subpictures  together  with  a 16-sketch  summary  for 
the  group.  Figure  5 shows  a section  of  the  data  preparation.  File  SUBPIX 
is  given  the  local  name  TAPEl  and  read  under  control  of  the  following  Fortran 
statements: 

READ(1,1) ( (NPIX(I,J) ,1=1,80) ,J=1,64) 

1 F0RMAT  (2020) 

Figure  6 shows  the  entire  array  of  subpictures.  The  pictures  were  pro- 
duced, six  to  a frame,  on  the  4060  with  a separate  program  designed  to  aid 
in  the  job  of  correcting  errors  in  the  data.  The  program,  named  DISPLAY, 
appears  in  Appendix  A. 
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Figure  3.  Data  Card  Representing  a Portion  of  the  Terrain  Model  MAP 


ROADS  AND  BUILDINGS 


Figure  4.  Typical  Sketch  Used  to  Formulate  Subpictures 


1^'"  1 

r‘ 1 
1 

4^ 

■tfS' 

liniiKfM 

i A 

ih.' 

J / 

:■ 

li  !- 

y 

H 

In  use,  the  terrain  model  is  accessed  by  way  of  a subroutine  named  TLU 
(lYDS,  JXDS,  JGRAY) . It  is  entered  with  rectangular  coordinates  in  integer 
yards  measured  from  the  upper  left-hand  corner  of  the  model  (across  is  lYDS 
and  down  is  JYDS) . In  the  first  stage  of  the  table-look-up,  it  is  determined 
in  which  element  of  MAP  the  desired  point  lies,  and  the  subpicture  designa- 
tion number  is  retrieved.  In  the  second  stage,  the  proper  point  within  a 
subpicture  is  resolved,  retrieved,  and  unpacked  into  a single  digit  for  return 
as  variable  JGRAY.  This  subroutine  appears  in  the  appendix  together  with 
the  rest  of  program  VTRAJ3. 


OPTICS  SIMULATION 

Another  major  block  in  the  simulation  is  optics,  wherein  a subroutine 
named  PHOTO  (X,  Y,  Z,  VX,  VZ)  provides  the  function  of  rendering  a picture 
of  the  terrain  model  (possibly  including  a target)  as  seen  frcan  a point  above 
the  model  specified  by  coordinates  (X,  Y,  Z) . The  optic  axis  of  the  view 
is  assumed  to  be  along  the  vehicle  velocity  vector.  The  velocity  vector 
is  specified  by  VX  and  VZ  which  are  normally  negative  quantities.  The  program 
does  not  allcw  for  a Y-component  of  the  velocity  vector  or  for  a roll  angle 
for  the  camera  away  from  the  upright  position.  Figure  7 shows  the  various 
coordinates. 


Figure  7.  Coordinates  Used  for  the  Optics  Simulation 
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The  output  of  subroutine  PHOTO  is  an  unformatted  write  to  tape.  The 
picture  format  used  by  the  simulation  programs  generally  is  a 128  by  128 
matrix  in  which  each  element  is  a single  octal  digit  representing  gray  levels 
in  the  same  way  as  for  the  terrain  model.  In  order  to  conserve  memory,  the 
octax  digits  are  packed  20  to  a word.  The  unformatted  write  is  from  a two- 
dimensional  array  named  IPAK  (7,128).  Figure  8 illustrates  how  the  array 
is  interpreted  as  a picture.  Also  shown  in  this  figure  is  an  array  of  four 
variables,  IHDR(4),  which  is  not  written  out  by  PHOTO  itself,  but  is  expected 
on  tape  by  the  various  other  programs  that  use  the  tape  as  input.  A value 
of  IHDR(1)=0  IS  used  to  indicate  an  end-of-file  condition.  Also  illustrated 
in  Figure  8 is  the  format  used  when  such  data  is  transcribed  to  35mm  film. 
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Figure  8.  128^  Pixel  Data  Format  With  the  Corresponding  35mm  Film  Format 

Internally,  subroutine  PHOTO  geometrically  projects  a grid  (shown  in 
Figure  7)  having  an  array  of  128  by  128  square  apertures  onto  the  surface 
of  the  terrain  model.  This  grid  simulates  a charge  transfer  device  (CTD) 
television  sensor.  The  fieid-of-view  angle  may  be  controlled  by  the  para- 
meter FOCAL,  which  is  the  distance  from  the  grid  to  the  lens.  FOCAL  is  in 
units  of  length  equal  to  a side  of  the  grid.  This  parameter  is  not  read 
in  but  is  set  in  the  program  by  a DATA  statement. 

The  four  corners  of  each  projected  aperture,  or  a pixel  in  the  simulated 
CTD  sensor,  are  located  as  integer  coordinate  points  on  the  terrain  model. 
These  points  will  generally  describe  a trapezoid  (Figure  9).  A subroutine 
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named  TRAP  then  uses  subroutine  TLU  to  obtain  all  of  the  terrain  model  elements 
that  comprise  the  trapezoid  and  forms  an  average  of  their  values.  This  average 
becomes  the  gray  level  value  of  the  corresponding  pixel  in  the  array  IP^. 

Since  the  sides  of  the  trapezoid  are  generally  not  parallel  to  the  y-axis 
of  the  coordinate  system,  a variable  slope  "staircase"  algorithm  is  used 
to  obtain  an  approximation. 


CAHT  OF 

TEDIIAIN  MODEL 

3333555565 

3335555555555555 


Figure  9.  Trapezoidal  Projection  of  a Pixel  on  the  Terrain  Model 


trajectory  simulation 

A simple  subroutine  named  FLIGHT  (X,  Z,  VX,  VZ,  DELT)  is  used  to  update 
the  coordinates  (X,  Z)  of  the  projectile  at  a time  interval  DELT.  Since 
the  trajectory  is  assumed  to  have  no  Y-velocity,  this  parameter  does  not 
appear.  X,  Z,  and  velocities  VX  and  VZ  are  updated.  One  hundred  iterations 
are  performed  for  each  entry  of  the  subroutine.  The  equations  used  are 


dt  » At/100 


dt 
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where 


VXi+l  = VXi  + AV^ 

Av^  = - (32.16  + Cd  Vz^)  dt 


At  = DELT  = time  interval  between  pictures 
Vx  = X - velocity 
Vz  = Z - velocity 

Cj3  = drag  coefficient  (set  to  0.04) 

Figure  10  shows  a typical  trajectory  and  view  angle. 
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Figure  10.  Typical  Trajectory  and  View  Angle 
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MOVING  TARGET 


A simple  provision  in  the  basic  terrain  model  subroutine  TLU  is  used 
to  superimpose  a black  square  moving  object  onto  the  scene.  The  size  and 
location  of  the  target  is  determined  by  the  three  variables  (IT,  JT,  and 
NSIZE)  which  are  in  CC»4M0N  storage.  NSIZE  is  the  size  of  the  square  target 
(given  in  yd^) . The  coordinates  IT  and  JT  are  updated  by  a subroutine  CRAWLR 
which  simply  updates  the  position  of  the  square  according  to  straight-line 
velocities  and  the  same  time  interval  used  by  subroutine  FLIGHT.  The  actual 
super  imposition  occurs  in  subroutine  TLU. 

A full  Fortran  listing  of  program  VTRAJ3  and  a set  of  operating  instruc- 
tions are  found  in  Appendix  A.  Figure  11  shows  a typical  video  data  file. 
There  are  39  pictures  taken  at  0.12-sec  intervals  along  the  trajectory  of 
Figure  10.  A column  of  header  information  is  beside  each  picture.  These 
variables  are  (starting  at  the  top)  Frame  Number,  X-coordinate  (ft),  Y- 
coordinate  (ft) , and  Z-coordinate  or  altitude  (ft) . A moving  target 
is  shown  in  the  pictures.  It  is  a black  square,  15  ft  on  a side,  moving 
up  and  to  the  right  at  a speed  of  29  mi/hr.  It  appears  in  the  upper  left- 
hand  corner  of  frame  33  and  can  be  traced  from  there. 

The  output  file  of  VTRAJ3  is  transcribed  into  Stromberg-Car Ison  4060 
(and  then  to  photographic  format)  by  a program  named  MOVIE.  The  listing 
and  instructions  for  this  program  also  appear  in  Appendix  A. 


TARGET  DISCRIMINATION 


Once  data  files  containing  sequences  of  pictures  (such  as  those  of  Figure  11) 
are  established,  the  next  task  is  to  find  a process  that  will  discriminate 
moving  objects  from  surrounding  stationary  clutter.  If  the  observation  of 
a target  area  could  be  made  from  a stationary  position,  moving  target  discrimina- 
tion would  be  very  easy.  It  would  be  necessary  only  to 

1.  Take  a picture 

2.  Store  it 

3.  Wait  a suitable  time  interval 

4.  Take  another  picture 

5.  Compare  the  new  picture  and  the  stored  picture,  pixel  by  pixel. 

Those  objects  in  motion  would  show  up  prominently  in  the  ccxiqparison  process. 
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Computer-Generated  Pictures  of  an  Area  Containi 
Moving  Target 


-Generated  Pictures  of  an  Area  Ci 
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Figure  11.  Computer-Generated  Pictures  of  an  Area  Conta 
Moving  Target  (Continued) 


Computer-Generated  Pictures  of  an  Area  Containing 
Moving  Target  (Continued) 


An  appreciation  for  the  nature  of  the  required  image  transformation  [ 

required  can  be  gained  from  Figure  7.  The  picture  taken  is  that  of  a trape-  1 

zoidai  "footprint"  on  the  ground.  As  the  trajectory  progresses,  motion  along  ' 

the  optic  axis  will  cause  the  footprint  to  shrink  to  a similar  but  smaller  ^ 

trapezoid.  The  effect  upon  the  picture  is  an  enlargement  and  an  apparent  | 

motion  away  from  the  optic  axis  of  all  objects  in  the  field  of  view.  Some  )| 

objects  will  disappear  from  view  and  some  fine  details  will  be  resolved  for  S 

the  first  time.  | 

There  are  other  components  to  the  observer's  motion.  While  rotation  I 

about  the  optic  axis  has  been  ruled  out,  the  optic  axis  can  still  pitch  over  t 

to  a steeper  angle  of  descent.  The  main  effect  from  a pitching  motion  is 
an  upward  vertical  displacement  of  all  the  objects  in  the  picture.  This 

shift  is  accompanied  by  second-order  distortions  that  result  from  the  fact  ; 


that  pitch  causes  the  new  trapezoidal  footprint  to  be  sanewhat  dissimilar 
from  the  original.  This  dissimilarity  means  that  the  enlargement  of  the 
scale  of  the  picture  is  not  uniform. 

For  simulation  purposes,  the  second-order  effects  are  ignored.  An  algorithm 
for  them  has  not  been  devised,  and  such  a step  may  prove  unnecessary.  The 
image  transformation  device  (Figure  12)  is  capable  of  enlarging  a picture 
and  shifting  it  upward.  Two  considerations  that  remain  are 

1.  How  are  standard  electronic  components  to  be  arranged  to  carry  out 
these  transformations? 

2.  How  are  the  control  parameters  (how  much  enlargement  and  how  much 
shift)  to  be  generated? 

First,  the  electronics  and  their  computer  simulation  will  be  treated. 

While  a microprocessor  is  assumed  to  be  a part  of  the  electronic  mix,  it 
cannot  be  reasonably  assumed  that  a microprocessor  (or  even  a fairly  large 
computer)  capable  of  handling  the  quantity  of  data  generated  by  a television 
sensor  will  appear  in  the  near  future.  The  microprocessor  in  the  system 
simulated  in  this  study  will  be  a "supervisor"  that  performs  numerical  com- 
putations in  support  of  other  electronic  devices  that  actually  process  video 
data. 


The  system  simulated  assumes  that  video  information  is  generated,  stored, 
and  processed  in  long  analog  shift  registers.  Given  such  a storage  medium 
for  video  data,  image  shifting,  horizontal,  and  vertical  is  easily  accomplished 
by  supplying  the  storage  shift  register  with  extra  shift  pulses  after  the 
entire  image  has  been  read  in.  An  image  magnification  process  will  prove 
to  be  more  trouble. 

Figure  13  illustrates  the  processing  needed  when  an  image  to  be  magnified 
is  stored  in  a long  shift  register.  The  upper  portion  of  the  figure  contains 
a 16  X 16  matrix  (folded  256-bit  shift  register)  containing  the  digitized 
image  of  the  letter  "A".  Superimposed  on  this  matrix  is  another,  slightly 
smaller,  16  x 16  matrix  representing  the  area  that  is  to  be  expanded  into 
the  full-size  picture  format.  The  lower  portion  of  Figure  13  shows  the  result 
of  the  magnification  process.  Each  pixel  in  the  output  matrix  is  computed 
from  the  weighted  contents  of  up  to  four  pixels  in  the  original  matrix. 

The  simulation  uses  a subroutine  with  the  Fortran  name  of  MAGFY  (IPAK, 
JPAK,  ISHIFT,  JSHIFT,  MAGGIE).  The  purpose  of  this  subroutine  is  to  simu- 
late the  image  transformation  device  (Figure  12) . Its  action  is  to  fetch 
a picture  from  IPAK  (7,128),  shift  it  in  the  I (horizontal)  direction  by 
ISHIFT  pixels,  shift  it  in  the  J (vertical)  direction  by  JSHIFT,  increase 
its  scale  by  an  amount  controlled  by  MAGGIE,  and  place  the  result  in  JPAK 
(7,128). 
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15^  Segment  16 ^ Segment 


Figure  13.  Sample  Magnification  Process  of  a 15^  Segment  Expanded 

into  a 16^  Segment 

Figure  14  illustrates  the  magnification  process  on  a pixel  level.  It 
shows  two  superimposed  grids  representing  the  pixel  array  of  an  input  picture 
(the  larger  and  darker  grid)  and  the  pixel  array  of  an  output  or  magnified 
grid.  The  magnification  occurs  as  the  contents  of  the  smaller  squares  are 
reinterpreted  in  a framework  that  has  the  same  dimensions  as  the  input  grid. 
The  process  by  which  the  contents  of  a given  output  pixel  is  formed  follows; 

I,J  = indices  specifying  the  pixel  being  filled  in  the  output  picture 

i,j  = indices  specifying  the  pixels  being  accessed  in  the  input  picture 

f(I,J)  = contents  of  pixel  (I,J) 

g(i,j)  = contents  of  pixel  (i,j) 

M = magnification  to  be  applied 

Then  to  compute  f(I,J)  (illustrated  by  the  shaded  square  in  Figure  14), 
the  machinery  simulated  or  the  computer  algorithm  must  find: 


i = integer  part 


j * integer  part 
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Figure  14.  Superimposed  Pixels  From  an  Input  Picture  and  a 
Magnified  Output  Picture 

M • ^fractional  part  of 
1 - T 

M • j^fractional  part  of 
1 - V 

J)  = U*W*g(i,j) 

+ T*Wg(i,j+l) 

+ U‘Vg(i+l,j) 


+ T*V*g(i+l, j+1) 


In  the  shaded  area  illustrated  in  Figure  14: 


(I,J)  = (4,5) 

M = 1.11111... 

(i.j)  = (3,4) 

T = 0.66666... 

U = 0.33333... 

V = 0.55555... 

W = 0.44444... 

Figure  15  outlines  a possible  configuration  of  analog  shift  (CCD)  registers 
and  high-speed  sequencing  logic  to  effect  the  above  process.  Figure  16  shows 
a sequence  of  pictures  that  have  been  processed  by  subroutine  MAGFY.  The 
picture  at  top  center  was  used  as  input.  The  sequence  of  pictures  below 
show  the  results  of  applying  various  control  parameters.  The  column  of  param- 
eters beside  each  of  the  transformed  pictures  represent  sequence  number, 
horizontal  shift,  vertical  shift,  and  magnification  (128/MAGGIE). 

It  can  be  seen  that  as  the  shift  quantities  are  held  constant,  the  in- 
creasing values  of  magnification  (decreasing  128/MAGGIE)  cause  the  scale 
of  the  picture  to  expand  and  the  features  to  enlarge.  This  sequence  bears 
some  outward  resemblance  to  earlier  sequences  wherein  each  new  picture  was 
rederived  from  the  terrain  model. 

With  simulations  in  hand  for  the  various  black  boxes  of  Figure  12,  the 
next  step  undertaken  was  to  develop  a process  for  supplying  the  proper  values 
of  magnification  and  shift  parameters  to  the  image  transformation  device. 

These  values,  which  must  properly  compensate  for  vehicle  motion,  should  not 
depend  upon  specific  internal  information  of  the  trajectory  being  followed. 

Such  data  could  only  be  supplied  by  a data  link  or  by  expensive  internal 
sensors . 

The  process  evolved  in  the  simulation  effort  is  one  that  finds  proper 
values  for  the  parameters  by  experimentation.  Two  quantities,  vertical  shift 
and  magnification,  are  optimized  for  the  best  vehicle  motion  compensation 
attainable.  The  quantity  used  to  score  the  various  experimental  results 
is  the  correlation  between  the  stored,  transformed  picture  and  a picture 
taken  after  a time  interval  (0.12  sec). 
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Figure  15.  Hardware  Array  to  Effect  Image  Transformation 
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Correlation  (j)  is  an  index  by  which  two  pictures  are  compared  pixel  by 
pixel.  Usually,  correlation  is  normalized  so  that  4>  = 1.0  for  two  perfectly 
matched  pictures;  = 0 is  an  improbable  result;  and  4>  = 0.5  is  the  expected 
value  for  two  unrelated  pictures.  If  the  correlation  of  two  similar  pictures 
as  a function  of  displacement  is  determined,  typically  one  finds  that  the 
correlation  function  has  the  shape  of  a cusp,  with  the  peak  occurring  at 
the  displacement  that  provides  the  best  overall  match.  A similar  function 
is  found  with  respect  to  scale  or  magnification.  The  exact  shape  of  the 
curve  IS,  of  course,  dependent  upon  the  contents  of  the  picture.  Some  kinds 
of  pictures  (e.g.,  ones  containing  a marked  spatial  periodicity)  can  greatly 
modify  the  correlation  function  from  the  typical  curve.  It  is  assumed  in 
this  simulation  effort  that  such  pictures  are  too  atypical  to  be  taken  into 
account . 
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Figure  16.  Typical  Picture  (Top)  Operated  Upon  by  Various  Values 

of  Magnification 


The  simulated  process  for  automatically  designating  targets  resides 
in  a program  with  the  FortrcUi  name  XMABK.  This  program  accepts  target  data 
files  from  program  VTRAJ3  and  experiments  with  displacements  and  magnifi- 
cations so  as  to  find  the  best  correlation  between  a stored,  transformed 
picture  and  a picture  taken  at  a later  time.  It  then  subtracts  the  best 
correlated  pair  of  pictures  it  can  find,  and  searches  the  difference  picture 
for  targets.  If  it  finds  a target,  it  superimposes  a pair  of  intersecting 
black  lines  over  the  target  in  the  input  picture  and  writes  it  out  as  an 
output  file.  The  program  also  provides  a printed  record  of  its  experimental 
progress  during  the  flight.  Figure  17  illustrates  the  overall  data  flow 
of  a simulated  trajectory  and  target  designation. 


Figure  17.  Data  Flow  of  Simulation 


The  actual  process  within  XMARK  is  a rather  crude  one  which  represents 
a beginning  for  an  evolutionary  process.  Actually,  scheduling  the  processing 
implied  in  real  hardware  might  prove  difficult.  It  is  also  expected  that 
due  to  further  refinements,  such  scheduling  will  also  prove  unnecessary. 

Figure  18  is  a flow  chart  of  XMARK.  Initially  the  program  reads  in  one  "old" 
picture  from  the  target  file  and  one  card  containing  a beginning  value  for 
the  two  critical  parameters  JJO  = vertical  shift  and  MAGO  = magnification. 

The  program  then  enters  into  a loop  wherein  it  reads  in  a "new"  picture  from 
the  target  and  enters  the  sequence  of  computations  in  Figure  18.  There  are 
calls  to  several  subroutines  whose  functions  are  outlined  below. 

Subroutine  MAGFY,  which  has  already  been  described,  simulates  the  image 
transformation  device.  Its  calling  sequence  contains  three  parameters; 

ISHIFT,  or  horizontal  shift  in  integer  pixels;  JSHIFT,  or  vertical  shift 
in  integer  pixels;  and  MAGGIE,  or  magnification  factor.  In  program  XMARK, 
ISHIFT  is  always  zero  since  yaw  in  the  projectile  has  been  disallowed. 

Subroutine  C0RR2  computes  the  correlation  between  two  128  by  128  pic- 
tures. It  adds  the  absolute  values  of  the  differences  between  the  contents 
of  corresponding  pixels  in  two  input  pictures.  It  returns  with  the  raw  summa- 
tion as  SUM,  and  also  ccanputes  a normalized  correlation  COEF  which  reaches 
1.0  for  two  identical  pictures.  Only  the  raw  summation  SUM  is  used  in  pro- 
gram XMARK. 

Subroutine  CURFIT  is  a curve-fitting  algorithm  for  use  with  trial  values 
of  correlation.  Since  correlation  as  a function  of  either  magnification 
or  shift  tends  to  be  a cusp,  this  subroutine  attempts  to  determine  the  best 
value  of  those  parameters  by  finding  the  peak  of  that  cusp.  Figure  19  shows 
a typical  correlation  cusp  as  a function  of  some  independent  variable.  In 
practice,  that  independent  variable  can  only  be  evaluated  for  integer  values 
given  the  quantized  nature  of  the  pictures. 

The  subroutine  receives  three  values  of  correlation  from  three  evenly- 
spaced  values  of  the  independent  variable  i.  It  computes  the  coefficients 
for  a second-order  polynomial  that  fits  those  three  points,  and  then  computes 
the  second  derivative  of  that  function  at  the  center  value  of  i.  If  the 
second  derivative  is  positive  (as  it  would  be  for  the  three  values  desig- 
nated by  X on  Figure  19) , it  computes  the  inter r t of  the  polynomial  and 
the  line  where  correlation  =1.0  and  returns  an  incremental  value  for  i equal 
to  Aj,  (as  shown).  If  the  second  derivative  is  negative  (as  it  would  be  for 
the  values  designated  by  "0"  on  Figure  19),  then  the  subroutine  finds  the 
nearest  integer  value  of  i for  which  the  parabolic  curve  has  a null  first 
derivative  and  again  returns  an  incremental  value  of  ^2  shown). 

Subroutine  MINUS  simply  subtracts  two  128  by  128  pictures  to  produce 
a difference  picture.  The  differences  of  the  values  of  the  pixel  pairs  are 
rescaled  to  have  the  range  (0,7). 
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Figure  19.  Correlation  Function  Curve  Fit 

Subroutine  SIFT  searches  through  a difference  picture  looking  for  a 
3x3  matrix  of  values  that  lie  outside  the  range  (3,4).  When  it  finds  one, 
it  generates  two  intersecting  black  lines  over  the  spot  and  superimposes 
i them  over  the  picture  that  is  the  minuend  of  the  subtraction  process  in  MINUS. 

' Subroutine  SIFT  also  reports  its  find  on  the  printer. 

' The  main  simulation  program  XMARK  uses  these  subroutines  (plus  two  other 

housekeeping  subroutines)  in  a sequence  (Figure  18) . The  program  first  reads 
in  one  picture  from  the  target  file  and  reads  a card  containing  starting 
values  for  vertical  shift  (JJO)  and  magnifications  (MAGO) . It  then  enters 
a loop  during  which  it  reads  in  a subsequent  picture,  updates  JJO  and  MAGO 
for  the  best  correlation  between  the  modified,  previous  picture  and  the  current 
picture,  subtracts  the  two,  and  searches  the  difference  for  targets. 

The  first  step  in  optimizing  the  values  of  JJO  and  MAGO  is  to  trans- 
form and  correlate  the  previous  picture  with  the  present  one  three  times, 
using  JJO-2,  JJO,  and  JJ04-2.  These  three  correlations  are  processed  by  sub- 
routine CURFIT  in  order  to  obtain  a corrected  value  of  JJO.  Using  this  cor- 
rected JJO,  the  program  then  transforms  and  correlates  three  more  times  using 
values  of  magnification  corresponding  to  no  magnification,  the  previous  magnifi- 
cation, and  twice  the  previous  magnification.  These  results  are  also  pro- 
cessed by  CURFIT  to  obtain  a possible  correction  to  MAGO.  If  the  correction 
obtained  is  nonzero,  another  transformation  is  done.  The  final  optimized, 
transformed,  previous  picture  is  subtracted  from  the  current  picture  and 
the  result  searched  for  targets  by  subroutine  SIFT.  During  this  procedure. 
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a printout  occurs  showing  the  values  of  JJO,  NAGO  used,  the  type  of  optimi- 
zation undertaken  (S,  M,  or  R for  shift,  magnification,  or  recomputation) , 
and  the  three  values  of  correlation  (displayed  as  the  raw  summation)  obtained. 
If  a target  is  found,  a print  line  identifying  the  coordinates  appears. 

Figure  20  shows  a sample  picture  output  from  program  XMARK  operating 
on  the  target  file  displayed  in  Figure  11.  In  the  25th  frame,  representing 
an  altitude  of  1789  ft,  the  simulation  program  correctly  identified  and  desig- 
nated the  object  in  the  field  of  view  that  was  in  motion. 


CONCLUSIONS 


Several  runs  of  the  simulation  were  made.  The  trajectory  was  changed 
only  in  landing  area;  the  size  of  the  target  was  maintained  at  15  ft^  and 
its  speed  was  maintained  at  30  mph.  The  target  direction  on  the  ground  was 
varied,  however.  Figure  21  shows  the  designations  made.  Several  false  desig- 
nations occurred.  These  were  always  prominent  objects  found  near  the  edge 
of  the  picture.  Assuming  that  designations  near  the  edge  of  the  picture 
are  to  be  ignored,  the  targets  were  properly  designated  at  the  following 
altitudes: 


Run 

Simulated 
Altitude  (ft) 

Extrapolated  for 
512  Lines-Altitude 

1 

1789 

7156 

2 

652 

2608 

3 

1030 

4120 

4 

904 

3616 

5 

1283 

5132 

Average 

1132 

4528 

The  pixel  array  size  simulated  was,  of  course,  only  128  x 128 — a small 
number  for  CCD  sensors.  If  the  array  size  were  increased  to  a state-of-the- 
art  figure  like  512  x 512,  the  altitudes  at  which  designations  could  be  ex- 
pected would  rise  proportionally  (by  similar  triangles)  to  values  four  times 
those  listed.  Similarly,  the  expected  altitude  for  designations  would  be 
directly  proportional  to  target  size  and  target  speed.  Generally,  targets 
in  the  upper  portion  of  the  picture  are  harder  to  detect  than  those  in  the 
lower  portion,  because  perspective  makes  them  appear  smaller. 
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In  the  course  of  these  simulations,  an  alternate  hardware  configuration 
has  been  conceived.  The  main  advantage  of  this  alternate  is  the  elimination 
of  the  hardware  for  image  transformation.  When  it  was  assumed  that  the  image 
transformation  device  would  be  used,  the  parameters  that  were  continuously 
optimized  during  the  flight  were  shift  and  magnification,  with  the  time  between 
picture  "takes"  (AT)  held  constant.  If  magnification  were  to  be  kept  at 
some  fixed  value,  presumably  one  could  optimize  shift  and  AT  instead.  With 
magnification  a constant,  it -would  be  possible  to  obtain  that  function  by 
optical  means.  Figure  22  shows  a hardware  reconfiguration  that  takes  advantage 
of  these  possibilities.  The  image  transformation  and  the  separate  image 
storage  blocks  are  eliminated.  The  separate  image  storage  is  eliminated 
because  the  second  image  sensor  which  is  added  to  accept  the  optically  mag- 
nified image  can  also  serve  as  a storage  and  shifting  device. 


t 


At 


Figure  22.  Alternate  Hardware  Configuration  to  Eliminate  Image 

Transformation  Device 


This  second  configuration  has  not  been  simulated  because  of  program 
deadlines;  however,  a future  simulation  is  anticipated. 
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APPENDIX  A 


PROGRAM  LISTINGS 

PROGRAM  VTRAJ3 

Program  VTRAJ3  creates  files  picture  data  given  terrain  model  data, 
trajectory  data,  and  target  data. 

Scope  and  Data  Cards  for  VTRAJ3 

ATTACH,  TAPEl,  SUBPIX,  ID..= 

ATTACH,  TAPE2,  MAP,  ID  = 

REQUEST,  TAPE4,  *PF. 

ATTACH,  VTRAJ3,  ID= 

VTRAJ3. 

CATALOG,  TAPE4,  XXX,  ID=... 


(■^89) 

TRAJECTORY  CARD 
TARGET  CARD 
(6789) 

Where  TRAJECTORY  CARD  contains  X,  Y,  Z,  VX,  VZ,  AT,  N in  Format 
(6F10.0,  14) 

X = Initial  x coordinate  of  projectile  in  feet 
Y = Initial  Y coordinate  of  projectile  in  feet 
Z = Initial  Z coordinate  of  projectile  in  feet 
VX  = Initial  X-velocity  in  feet/sec  (negative  number) 

VZ  = Initial  Z-veiocity  in  feet/sec  (negative  number) 


N = number  of  pictures  to  be  taken 

and  where  TARGET  CARD  contains  TS,  TY,  TVX,  TVY,  SIZE  in  Format 
(5F10.0) 

TX  = Initial  target  X-position  in  feet 

TY  = Initial  target  Y-position  in  feet 

RVX  = Target  X-velocity  in  feet/sec 
TVY  = Target  Y-velocity  in  feet/sec 
SIZE  = Target  size  in  yards  square 

Subroutine  Crawler — Updates  the  position  of  the  target  on  the  terrain 

map  during  flight 

Subroutine  Photo — Performs  iteration  at  128  x 128  pixels;  projects  each  pixel 

on  terrain  model  and  finds  corners  of  each  resulting  trapezoid 

Subroutine  Flight — Updates  the  position  of  the  projectile  during  the  flight 

Subroutine  TLU — Does  table- look-up  on  terrain  model  given  coordinates  in 
yards 

Subroutine  Trap — Accepts  location  of  four  corners  of  a projected  pixel  on 

the  terrain  model  and  computes  average  gray  level  for  pixel 

Subroutine  Pakrat — Accepts  an  array  containing  one  raster  line  in  128  locations 
and  packs  the  information  into  seven  locations  within 
a 7 X 128  array 


PROGRAM  MOVIE 

Program  MOVIE  accepts  files  of  picture  information  and  renders  them  on 
35mm  film  using  the  Stromberg-Car Ison  4060. 

Scope  cards  for  MOVIE 

ATTACH, MOVIE, ID=. . . . 

ATTACH, SYSLIB. 

LIBRARY, SYSLIB. 

ATTACH , TAPE4 , XXX , ID= 

REQUEST,  TAPE3*PF. 
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MOVIE. 


DISPOSE,  TAPE3,  PR=CCB. 

(Sg) 

PROGRAM  XMARK 

Program  XMARK  simulates  target  designation  hardware  operating  on  simulated 
video  data  from  program  VTRAJ3. 

Scope  and  Data  cards  for  XMARK 


ATTACH, XMARK, ID = 


ATTACH, TAPEl, XXX, ID= (from  VTRAJ3) 

REQUEST , TAPE4 , *PF . 


XMARK. 


(TAPE4  now  contains  a copy  of  the  picture  file  on  TAPEl  with  the  possible 
super  imposition  of  cross-hairs  cn  designated  targets;  it  may  be  viewed  by 
means  of  program  MOVIE  if  desired.) 


(■^89) 

DATA 

(6789) 

Where  Data  Card  contains  n,  JJO,  MAGO  in  * format 
N = number  of  frames  to  be  processed 
JJO  = initial  guess  of  vertical  shift  (1) 

MAGO  = initial  guess  of  magnification  (126) 

Subroutine  MINUS — Subtracts  two  picture  arrays  pixel  by  pixel  and  places 
rescaled  difference  into  a third  array 

Subroutine  MAGFY — Receives  a picture  in  one  array,  shifts  it  by  specified 

numbers  of  pixels  in  I,J  directions,  applies  the  magnifi- 
cation specified,  and  places  the  result  in  a second  array 

Subroutine  CURFIT — Receives  three  values  of  correlation  (as  raw  scores)  that 
are  assumed  to  result  from  three  values  of  some  parameter 


spaced  apart  by  some  delta;  subroutine  canputes  correction 
increment  to  be  applied  to  parameter 

Subroutine  SIFT — Searches  through  a difference  pict  e array  for  a 3 x 3 

area  containing  gray  levels  outside  the  range  (3,4);  search 
is  left  to  right  and  down;  the  first  one  found  is  designated 
by  a cross  placed  over  the  picture  in  a second  array 

Subroutine  FETCH — Obtains  a single  unpacked  pixel  value  from  a picture  array 

Subroutine  C0RR2 — Computes  the  correlation  between  two  pictures;  result  is 
presented  as  a raw  score  as  well  as  a normalized  value 


PROGRAM  DISPLAY 

Program  DISPLAY  renders  the  file  of  subpictures  onto  film  via  the  Stromberg- 
Carlson  4060. 


PROGRAM  VTRAJ3 


3?0SRAM  /TRAJ3  7^/7-  0PT=1 


1 


5 


1 3 


15 


20 


25 


30 


35 


PROGRAM  \/Tk,AJ3  (INPUT, output, TAPei,TAPE2,TAP54) 
COMMON  NPIX (6 0,64 », IMAP (48,120 ), IT, JT, NS IZE 
01MENSI3S  i,H0R(4> 

1C0UNT='I 

PEAO  (1,1)  ((NPlK(i,JJ,I=l,eC» ,J=1,64) 

1 FORMAT  (2320) 

READ  (2,2)  ((IMAF(I,J) ,1=1,24) ,J=1, 120) 

READ  (2,2)  ((IMA3(I,J),I=25,48),J=1,120) 

2 FORMAT  (2402) 

READ  3,X,f ,Z, VX,tfZ,OELT,N 

3 FORMAT  ( 6Fi0, 0,14) 

C • • • 

3 ENTER  TARGET  POSITION  X ,V  FT , /^ELO ( VX  ,V Y ) /Si ZE 

C • • • 

READ  3,TJ,T1,)/J,/I,SIZE 

NSIZ£=SI2E-1. 

iT=TI/3. 

JI=TJ/3. 

4 ICOUNT  = rCOUNT4-l 
IHOR(l)=ICOUNT 
IN0R(2)=X 
IH0R(3)=Y 
IHDR(4)=Z 
WRITE  (4)  IHOR 

PRINT  9q,S,X, Y,Z,VX,VZ 
9?  FORMAT  (I7,3X,5F12.2) 

CALL  PHOTO  (X,Y,Z,VX,VZ) 

N = N-1 

IF(N.LE. 0) GO  TO  5 

CALL  FLIGHT  ( X , Z , VX , V Z , 0£ L T ) 

IF  (Z.LT, 120.0)  GO  TO  5 
CALL  CRAWLR  ( T1 , T J , tf I , V J , OELT ) 

GO  TO  4 

; IHOR(l)=Q 

WRITE  (4)  IHOR 

STOP 

ENO 


suaouriNE  :pahlr  74/74  opt=i 


1 SUBROUTINE  drawer  ( TI , T J , V I, V J , DEL T ) 

COMMON  NPIX(8r,64) ,IMFA(4e,l20) ,IT,JT,NSIZE 

c , , , 

C UPOATF.  TARGET  POSITION 

5 C ... 

TI=TI4VI*DELT 
TJ=Tj4\/J*DfLT 
IT  = TI/3. 

JT=TJ/3. 

10  RETURN 

END 
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SUBiDUTINE  3iOTO  7it/74  0PT  = 1 


t 


\ 

I 


1 


5 


13 


15 


20 


25 


30 


35 


40 


45 


E 0 


E 5 


SUBROUTINE  PHOTO  « X , V , Z , V X , VZ ) 

DIMENSION  10(129) , IN ( 1 29 ) . IPR ( 12 8 ) 

1 ,IPAK(7, 12P) 

DATA  F0CAL/5.C/.HAFPI/1. 5707963/ 

V • • • 

C COMPUTE  BOkESIGHT  ANGLE  THETA 

IFUZ.NE.C.O)  THETAsATANI  VX/WZ) 

IF  (y/Z.EQ.O  . 0)  THETAsHAFPI 

^ • • • 

: COMPUTE  initial  etao.phio  based  on  0EL=.5 

c 

eTAO=ATAN( .5/FOCAL) 

phio=theta+etao 

c • • • 

C COMPUTE  INITIAL  X-INTEP-CEPT  , J1  POINTER 

C • • • 

XIO=X-Z*TA'.(PHIO) 

Jl=XlO/3. 

if(xio.lt.:.o) ji=-i 

IF(PHIO. GE. HAFPI) Jl=-2 

C COMPUTE  INITIAL  ARRAY  OF  V-InTERCEPTS  I0(I) 

c • • • 

00  200  1=1.129 

RI=I 

EPS=-.5» (Ri/128. ) 

ylO=tPS*Z*:OS(ETAO)/(FOCAL*C3S(PHlO)) ♦Y 
2?:  IO(I)=YIO/3. 

C DOUBLE  DO,  COMPUTE  128  COLS  PEP  128  LINES 

Z ... 

00  30  0 J=l,128 

RJ=  J 

0ELS.5-R J/128. 

ETAsATAN (DEL/FOCAL) 

PHl=TH£TAf • TA 
XI*X-Z*T AN (PHI) 

J2=X  1/3. 

IF(XI.LT.O.O) J2=-l 
IF<PHI.GE. HAFPI) J2=-2 
00  400  1=1,129 
RI=I 

£PS=-.5^(RI/128.) 

YIN-EPS»Z*COS (ETA)/)FOCAC*COS(PHI) ) ♦Y 
43?  IN<I)=VIN/3. 

DO  500  1=1,128 
111=10  (I) 

112  = 10  (1*1) 

I21=1N(I» 

I22  = IN  (1*1) 

CALL  TRAP  < Jl, 32,111,112, 121,122, IG) 

5JD  IPR(I)=IG 

w ... 

C SAVE  scan  line  WITHIN  IPAK 

C ... 

call  PA<RAT  (J,1PR,IPA<) 
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PROGRAM  VTRAJ3 


SJB?3UTINE  ^HOTO 


0PT  = 1 


60 


c 5 


T 0 


C • • • 

: J2  SECOIFS  Jl»  INIU  BECOMES  IO(I> 

c • • • 

Jl=  J2 

00  600  1=1,129 
6j0  I0(1I=IN(I) 

c • • • 

C LOOP  ON  SCAN  LINES 

C • • • 

3CJ  CONTINUE 

C • • • 

C WkITE  OUT  PICTURE  ON  TAPE 

WRITE  (4)  IPAK 

RETURN 

END 


SUBiOUTINE  FLIGHT  74/74  OPT=l 


1 SUBROUTINE  FLIGHT  ( X , Z , X, tfZ, OELT) 

DATA  OT/.Ol/,  DRAG/. 04/ 

T = 0, 

1 0EL^X=-rfX*09AG*0T 

5 X=X»iVX*0ELVX/2. ) *DT 

VX=\/X*0EH/X 

OELVZ=-( 32.16+VZ*DRAG)*0T 
Z»Zt(VZ*0EH/Z/2.  »*DT 
VZ=VZ*OEH/Z 

13  T=T»0T 

IF(T.LT.OELT)GO  TO  1 

RETURN 

ENO 


SrMBOLIC  reference  HAP  (P=l) 

ENTRY  POINTS 
3 fli;tt 


l/ARI  ABLES 

SN  TYPE 

RELOCATION 

0 

DELT 

REAL 

F.P. 

35 

DELtfX 

real 

36 

OEL/Z 

peal 

31 

DRAG 

PEAL 

30 

OT 

peal 

34 

T 

PEAL 

0 

YX 

REAL 

F.P. 

0 

VZ 

real 

0 

X 

REAL 

F.P. 

0 

Z 

PEAL 

1 
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74/74 


0PT=1 


10 


15 


20 


25 


t 30 


35 


SUBROUTINE  TLU  ( I YOS , JT D S, JG RA Y) 

COMMON  SPIX(80,64),IMAP(4etl20) tlTiJT.NSIZE 

C TEST  IE  OM  TABULATED  REGION 

C • • • 

IF(IYOS.LT.O.OR.IYOS.GT.1S19)GO  TO  100 
IF( JYOS.LT. 0. OP. JYDS. GT. 4799)G0  TO  100 

A 

U . • • 

C TEST  IF  ON  target 

C ... 

IF(IYOS. GE.IT.ANO.IYOS.LE.  1 1 T ♦NS  I Z E» . AND . 

1 JYDS. GE.  JT. AND.  JYOS.LE.  (JT^NSIZEDGO  TO  ICl 

o ... 

C FINO  I,J,K  LOCAL  COOKOINAT^S 

C ... 

IM=  IYDS/40 
JH= JYDS/4P 

K=IMAP<I‘1^l,  JMfl)  ♦! 

I=IYDS-IM*- 0^1 
J=  JYOS-J'1*40  + 1 

C ... 

C LOOK  UP  ON  SUBPIX 

^ . 

LR=I/21 

NSHIFT=(I- (LR*20) )*3 
JR0H= ( (J-1) ’Z) ♦l^LR 
JPIX=NPIX( JPOM.K) 

JGRAV=S4IFT( JPIX.NSHIFT) ,AN0.7B 
RETURN 
C . 

C TREAT  OFF-MAF  SITUATION 

C • • • 

i:C  JGRAY=4 

IF( JYOS.LT. 2)  J3PAY=7 
RETURN 

1.1  JCRAY=0 
RETURN 
END 


1 1 


SYH30LIC  REFERENCE  MAP  (P=l) 


cNTLY  POINIi 
3 TLU 


variables 

51  I 

Sn  type 

INTEGER 

RELOCATION 

5b 

IM 

INTEGER 

1230  0 

IMA^ 

INTEGER 

ARRAY  / / 

25200 

IT 

INTEGER 

0 

I YDS 

INTEGER 

F.P. 

62 

J 

INTEGER 

3 

JGR4Y 

INTtG>^R 

F.F. 

57 

JM 

INTEGER 

56 

JPIX 

INTEGER 

65 

JROW 

INTEGER 

212)  1 

JT 

INTEGER 

/ / 

0 

JYDS 

INTEGER 

53 

< 

INTEGER 

63 

LR 

INTEGER 

0 

NPK 

INTEGER 

ARRAY  / / 

64 

NSHIFT 

integer 
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PROGRAM  VTRAJ3 


SL»3<JjriNE  ti;<ap  7‘*/7i,  nPT  = l 


1 SUBROUTINE  TRAP  ( J1 , J2 , 1 11 , II 2 , I 21 , 1 22 , TG ) 

JO= J2-J1*1 
JOP  = JO  ♦! 

IF(I21.Gfc. ill)  ISTAGl=JOP/ (I21-I11+1) 

5 IF«I2l.LT. Ill)  ISTAG1= JOP/ (I21-I 11-1) 

IF<  I22.GE. 112)  I STAGE = JOP/ (122 -1 12+1) 

IF (I22.lt, 112  I ISTAG2= JOP/ (122-112-1) 

IG=  Q 
N=C 

1C  DO  1 JO<:l,JO 

JYOS= Jl» JOX-l 
IL0LIM=I11 ♦(JOX+l)/ISTAGl 
IHILIM  = I12«-  ( JOX+1  )/ISTAG2 
00  1 lYOS  =ILOLIm,IHILIH 

15  CALL  TLU  (lYCJS,  JYDS.LEVLG) 

N=N  + 1 

1 IG=IG+LiKLG 

IG= ( CN/2) tlG) /N 
RETURN 

29  END 


SUBRDUTINE  PAKFAT  T4/7u  0PT=1 


1 SUBROUTINE  PAKRAT  (J,1PR,IPAK) 

OINENSION  IPR (126) ,IPAK(7,128) 

c • • • 

C CLEAR  OUT  J-TH  ROH  IPAK 

5 C 

00  1 K=l,7 

1 IPAK(K,J)  = '; 

C PACK  IFR(128)  INTO  IPAK  (K,J),K=1,7 

10  C • • • 

DO  2 N=l,12a 
HN= (N-1) /2C 
HNI  = MN+1 

ISHIFT=23- (N-EO’MN) 

15  KSHIFT=3*ISH1FT 

2 IPAK(NNI, J) =IPAK(MNI, J) , OR. SHIFT  ( IPR (N ) , KSHIF T ) 

RETURN 

END 


PROGRAM  MOVIE 


<»«?OGRAM  Hovre 


7i*/7u  0PT  = 1 


1 


7 


I'J 


1 j 


2J 


2> 


33 


35 


43 


•45 


63 


55 


pKOGkAh  movie  ( input, output ,TflP£3,TAPE4,TAPE51=0UTPUT) 
II  Me  NS  ION  IH0R(4t  , IPAK(7,12a»  ,LTAe(8)  ,I  PR  < 128  ) , 7 C ?0  C ) 
OATA  LTAB/lR4,lRP,lRC,lR3,lRQ,iR2,lRl,lR  / 

0 • • « 

C INITIALIZE  4063 

C • • • 

CALL  MODc:SG<Z,3) 

CALL  CRTIO  (7,5HNEZRH) 

CALL  PAGEG(Z,0,0,1> 

C StT  SCALE  TO  RASTEP  SIZE 

C • • « 

CALL  SETSMG(Z, 19,4095. > 

CALL  SETSMG(7,2J,3071. » 

C , , • 

C 4tA9  HI  ACER  - POSSIBLE  ENO  OP  RUN 

C • • • 

ii,  REA0(4)  IHOP 

IF(IHPP(H.NE,0)GO  TO  1 
CALL  PAGEG(Z, 0,0,1) 

CALL  CRTID  (Z,5H  ) 

call  EVITCMZ) 

STOP 

c • . • 

C WRITE  HEADER  INFORMATION  LARGE  CH.-UPPER  CASE 

C . . • 

1 CALL  SETSMG  (7,45,1.5) 

call  SETSMG  (Z, 52, 0.0) 

X=3i»i0. 

T=  2 JcO. 

CALL  FmTSG(Z,1,5,  J.IHORd)  , IPL) 

CALL  LEGN0G(Z,X,Y,5,IPL) 

Y=V-333. 

CALL  FMTSG(Z,1.5, 0,IHDP(2) ,IPL) 

CALL  L'  GNDGC Z.X  ,Y  ,5,IPL) 

Y=Y-333. 

CALL  FMTSG(Z,1,5,0,IHDR(3),IPL) 

CALL  L£GNOG(7,X,Y,5,IPL> 

Y=Y-333. 

CALL  FMTSG(Z,1,5,0,IHOP(«) ,IPL) 

CALL  LEGN0G(Z,X,Y,5,IPL ) 

c . . , 

C REAO  IN  COMPACTED  VIDEO  ARRAY,  INITIALIZE  PICTURE 

C . . • 

REAO  (4)  IPAK 

CALL  SFT5MG(Z,45,. 75 ) 

CALL  SETSMG(7,52, 2,0) 

X = 64, 

Y=3 JuA. 

C ... 

C LOOP  OVER  128  RONS 

C . . , 

00  111  J«l,128 

C ... 

C DECOMPOSE  J-TH  LINE  INTO  PRINT  ARRAY 

C ... 

00  1111  L*1,13 
1111  IPK  (L)*0 


44 
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PROGRAM  MOVIE 


PROGRAM  MOVie  74/7m  0PT=1 


6j 


6i 


7D 


75 


Of)  11  Nrl,12R 
'1N=  (N-1  J/2U 
)'N1  = MN  + 1 

^SHIFTS (N-20*MN)*3 

MYW  3k0=  (SHIFT  (IPAK(MNI  , J)  .NSHIFT  ) .and.  73)  +1 
ICH AK=LTA3 (MVMORO) 

M0X=  (N-1) /lO 
HOX I=Nnx»l 

KROT=6*(10-(N-10*NDX) ) 

II  IPF,  (NDXI)=  IPR(N0XI».0R,SHIFT(ICHAR,  KROT) 

c • • • 

C HRIT^  ON  CRT 

C « » a 

CALL  LCGNDG(7, X,Yal28, IPRI 

III  T=T-23, 

CALL  PAGEG(Z,0.0,H 

GO  TO  10 

£NU 


SYMBOLIC  1 

Rc.  FERt  NCc 

MAP  (R=ll 

ENTRY 

POINTS 

10  22  3 

MOVIE 

VARIABLES  SN 

TYPE 

RELOCATION 

US  6 1 

ICHAR 

INTEGEP 

10  5 65 

IHOR 

INTEGER 

lo57l 

IRAK 

INTEGER 

ARRAY 

10  5*1 

I PL 

INTEGER 

1240  1 

IPR 

INTEGER 

ARRAY 

10552 

J 

INTEGER 

IjS  64 

KR3T 

INTEGER 

10  5 53 

L 

INTEGEP 

12371 

LTAB 

INTEGER 

ARRAY 

10555 

MN 

INTEGER 

1J556 

MNl 

INTEGrR 

10560 

MYH0R3 

INTEGER 

10564 

N 

INTEGER 

10562 

NDX 

INTEGER 

10563 

NOXI 

INTEGER 

105  57 

NSHIFT 

INTEGER 

10547 

X 

REAL 

105  50 

Y 

real 

126U  1 

z 

REAL 

ARRAY 

FILE  NAMES 

MODE 

0 

INPUT 

2343  OUTPUT 

4106 

TAPE3 

2043 

TAPtSl 

EXTERNALS 

TYPE 

ARCS 

CRTIO 

2 

EXITG 

FMTSG 

6 

LEGNOG 

M03ESG 

2 

PAGEG 

SETSMG 

3 

INLINE 

FUNCTIONS 

TYPE 

ARCS 

SHIFT 

NO  TYPE 

2 INTRIN 

STATtMENT  LABELS 

U2iu  1 10236  10 

0 111  0 1111 
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ARRAY 


J 


■1  ..-rg 


^ tfN  J* 


OOO  OOO  iiOO  OOO  OOO  \0  ooo  c?cic>o 


PROGRAM 


IH.-  f «G£  iS-^EST  QUALITY  PBACIi..---4 
XMAI^jOiR  (JOr  X it UiU»lSHED  TO  DDC  ^ 


>?33RAN  XiiARK 


73/7^  0PT=1 


1 


PROGRAM  X>iAR<(  IWPUT,  OU  T»  UT  , ThPc  1 .TAPE 4) 

• • • 

PROGRAM  TO  OFTOmZL  COfRcLATION  AMO  WRiTE  OUT 
DIPPLRtNCE  PICIURE 


• • • 

0INEMSI04  IH0R<4)  .EPIX0(7,12b»  , NPI XT ( 7 . 128 ) . 
lNPiXN(7, 12o) ,NPIX0(7,  120) 

• • • 

REmO  in  CAt^Q  PARAMETERS  AnD  PRINT 

• • • 

READ  ♦ , N. JJO.MASO 
PRINT  9.N.JJ0,MA&0 

FOkMAT  (4H1  N=  .I4,>.X,hHJJ0=«I<*.4X,‘;hMAO0:.I4) 

PRINT  10 

j FORMAT (55rtt  I JJO  MAGO  X SUM-N  SUM-0 

PRINT  11 

1 FORMAT  (55M  - - - — 

2 FOkMAT  ( 1 5,  16 , 1 6 . 4X  , A1 , 1 9 , 1 11, 1 1 1 ) 


• • • 

READ  IN  OLj  PICTURE 

READ  (1)  IHDR 
IF  IIHDRCII.EQ. J)  STOP 
READ  (1)  MPIXO 

Initiate  outer  loop  computations 

00  7 1=1, M 

« « • 

REAO  IN  NEW  PICTURE 

• • • 

PEAO  (1)  IHOR 
IF  (IHCRdI.EQ.i)  STOP 
REAO  (II  NPIXN 


00  SET  OF  Three  shifts  and  correlations 

• • • 

JJ= JJO>2 

call  MAOFr  ( NPIXO.NPIXT.a. JJ.MAGO) 

C4LL  CORR2  ( NPIXN,NPIxr,COLF,SUMN) 
ISUHNsSJMN 

JJ> JJO»2 

CALL  MAGPV  ( NPIXOfNFIXT, J. JJ.MAGO) 

CALL  COkfi2  (NPIxn.nPIXT.ColF.SUMP* 

ISUMP=SUMP 

JJ=  JJO 

CALL  HAjFY  ( NPlX  0,  tIP  1 XT  . 0 , JJ.  M A GO 
call  CORR2  (NPIXN.NPaXT.COlF.SUMO) 
ISUM0=SJM0 

• • • 

PkINT  INTE'pmEOTAIE  result^ 

• • • 

X=1HS 

PRINT  12 .1 . JJU.MAOo.X  , ISUMN,I:>JM0.  ISUMP 

• • a 

DO  PARASOLaC  CURtlE  FIT 
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XMAKK 


PROGRAM  XMARK 


73/74  Oi’r:! 


IHIS  PAGE  IS  BEST  QUALi:  i PRACIICABLE 
FROW  COPY  RTRNISHE3  T'>  DOC  


CALL  CU^FII  ( SUIN, bUNO, S JMF, L.  0 t INCREJ 
JJO=JJOtI'JCRE 

no’sci  3F  rHi?EE  MAGNIFICATIONS  ANO  CO  i<RELA  T I ON  S 

• • • 

MA&=128 

call  MASFr  ( NPl/0, NPIXT , Uf JJO,  MAG) 

call  C0?F2  ( NPIXN.NPIX T.CQt F, iUMlO 

ISUMN=SUMn 

IZOOH=123-MAGO 

M«G=12fl-2*IZ00M 

call  MAGFY  (NPIXU.NPIXT.Q, JJ0,MAG) 
call  COiXRZ  ( NPIXN.nPIXT  ,C0LF,  SUMP) 

ISUMFsSUMP 
MAl.  = MAGJ 

CALL  MAGFY  ( NP I X 0, NPX X T , u , J JO,  M AG ) 

CALL  C0<?2  ( NPIXN,NPlXT,O0tF,SUM0) 

ISUMO= SUMO 

• • • 

PklNl  INTErtMEOIATL  RESULTS 

• • • 

X=1HM 

PRINT  12  ,I  , JJo.MAoO.X  , ISUMN,I:>UN0,  ISUMF 

t • • 

00  PARA33LiC  CURYE  FIT 

• • • 

OtLrA=12?,-MAGO 

CALL  CURFIT  ( SUMN , j JM 0 , S UMF , CEL T A , INCRE ) 

MAuO=MAGJ-INCRF 

IF ( HAGO. ST . 127)  MAG0=127 

• • • 

SKIP  RESJMP  IF  MAGNIFIlATION  NuT  CHANGES 

• • • 

IF  (INCRE. EQ.O)  SO  TO  A 

• • • 

RECOMPUTE  TRANSFOPMEO  PICTURE 

... 

CAlL  MhSFY  ( NPIX  0,  i.PI  XT,  Q , JJC,  MAGO) 

... 

PRINT  TRAC- 

• « • 

X=1HR 

PRINT  12,1, JJO, MAGO, X 

... 

SUjTkmCT  A'.RAYS  3Y  PIXEw  NFlxI -NFIX0  = NFIX0 

... 

CAlL  minus  (NPIx  T ,1,01  X b,  NPIxC) 

... 

SE«RCn  JUT  TARotT  ll<  SUSTRAlTED  ARRAY 

... 

CALL  ;>IFI  (NPI^kH,  NFIXD 

... 

M^ITc.  NEW  ►ICTURES  OUT  TO  HlE 


iM  uR ( 1 ) = I 
1M0R(2)  = J JS 
IH'JR(3)  s NASO 


I 

»J 


PROGRAM  XMARK 


'(MAiyK  73/7*.  0>T  = 1 


115 


120 


125 


13C 


rtSIIF  (4)  iHOR 
(h)  .mPIXN 

c • • • 

: TPAnbFC<  ^^FIXN  TO  nPIXC 

c • • • 

30  17  IDX=i,7 

00  17  JDX=1,12F 

17  NPI XC I lax, jOX ) = NPIXN  (lOx.JDX) 

c « • • 

; ENj  of  3UT!lK  loop  computatios 

c • • • 

7 CONTINUE 

w • • • 

C ARiTE  NJlL  NtAOt? 

U • • 9 

IHOR  <11  = 0 
WR1TE(4»  IHOP 
REMIND  4 
STOP 
ENO 


SUo7DJTI.<t  •■llNUS 


73/74  0PT=1 


1 


5 


10 


15 


20 


iUoROUTlNE  MINUS  ( NP I X T , NPi XC, N» I X D » 
jlMtNSION  NPI XT  <7,1^81 ,NPIX0 (7,1281 ,NPIX 0(7, 12d» 

C • • • 

C CLEAR  OUT  NPIXD  ARRAY 

• 

OG  2 1=1,7 
DO  2 J=i,128 
2 NPIXC(I,J>=0 

C , • • 

C subtract  TkO  PICTURES  AND  kESCALE  RESULT 

V • , , 

DO  1 1=1,7 
00  1 J=l,123 
00  1 K=4,E1,3 

NPIXLT=SHIFI (NPIXTtI , J) , K-1) . AND.7e 
NPiXL0=SMIFT(NPIX0(I,J),<-l).AN0.79 
NPIXLD=( 7»NPIXLT-NPIXlO) /2 
1 NPlxOd,  J»  =NPIXO<  I,  J)  .OR.SHiFr  (NFrXLD,bl“K» 

RETURN 
END 
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PROGRAM  XMARK 


SJa^JJTI'lE  'laGFY 


73/71-  0 = T=1 


1 J 


15 


25 


33 


35 


43 


45 
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55 


SUBkCUTINl  MiL.FY(IpHK,JPAK,i5nii^T,JSHlFT,-<AGGlL» 
JISENSIUlx  I-fiK(7,l^d)  ,JPaK(7.128I,IASSY(12«» 


bUuROUTlNc  Tc  SHIFT  FICIJkl  IN  IFAk  5Y  ( i SH IF T , Jb h I F T ) 
PIXtLS,  bNLARGt  ^LbULT  SY  - •.«  CN  I F I c»- T IG  H FAoTC?  G^ 
(1<;3/MAG),  AND  PlHLL  KlSULT  IN  JFAK. 


PPfcPROCSSS  HAG 

• • • 


MAG  = 126  - 2 • (126  - HAGuIL) 
M«GO=  («AG*lf’0C  : C)  /12l 
ImOPHF=123l  - JOO/MAG 
IShoR£  = ( 12{--M  AGI  /2 


LOOP  129  LxNES  ON  INDEX  J 


• • • 

DO  1111  J=ltl2b 


CLCAK  OJT  AiSEMGLY  ARRAY 


00  i 1=1, 12e 

1ASSY( II =C 


COMPUTt  J-JtFtN3ANT  OUANTITIES 


JN=  J*MAGO 
Jf,IP=JN/lLCOCO 
JNFP=JN-(  JNlP’UrC.JI 
Ji,FPNK=(  JNF=  *NGPMF)/lLPC3r 
JNFPS3  = 1 )0  f jn- JNPF.NM 


INNER  LJJF  128  PIXi-.LS  gN  iNGt  X I 


DO  111  1=1.128 


guMPUTc  AiijMlINj  FACTORS  P.J.R.S 


• • • 

IN=I*MA3D 

IMP=IN/  1C  I.C3  0 

INFP=1N- ( I NlP’luutCt ) 

INFPNMst INF  P*NCkMF I / IgC  0 Uu 

iNFFS6  = l-L  . i;0-lN  = Pl,M 

NP= ( INFPSa’JNFFSB) /IC  roc ) 

NQ= ( iNF3NP*  JNFPS9) / ICCOU  Q 

NR= (INF3S3*JNFFNMI/1CCuoC 

NS=  (INFPNM*  JNFPNM I /IL  JGg  J 


CQhPUTt  FETCH  INDICES 


ID  1= IN  IP  +1 SHPRl  + IShIF  T 
102=101+1 

jGl  = JNiP+I  bHF K”” ♦ J ShI F T 
JJ2=  JC  1*1 


call  S09»JJT1Nl  to  FlTLH  FIXl^S  FnOF  IfA< 


gAwg  FFTCh  ( IPAR,  IDl,  JDl  , K-*) 
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K'aJ'Jl  i -fUjudiHiij  XO  £OG 


PROGRAM  XMARK 


SJoiJUTlNt  -IfiOFY 


7J/7-  OPr=l 


uAL4.  FETGH  (IPAK,Iu2,  JUl.Ku) 
call  fetch  ( IPfiK, 101, JD2 , KP) 
call  fetch  < IRAK, IC2 , JD2 , KC) 

• • • 

COhPUTE  pixel  YA.Ut 

• • • 

L=  tKP*N=HKQ*NU  + Ki<*N<*KS*HS+5L  Jr-O  )/lCuaOO 
IF(l.GE.7)L=7 

• • • 

Place  in  AjSem3lv  aieray 

• • • 

lAiSY ( I»  =L 

• • • 

PACK  ASSPHiuY  ARRAY  INTO  JPA< 


CALL  PA<kAT(  J,lAbSY,  JPAO 

CONTINUE 

RETURN 

END 


SJaRDJTINE  PAKRAT 


73//-  0PT=1 


SUJROUTI„£  PAKRAT  tJ.IPR.IPAK) 

UIHENSION  iPP  ( l2o)  ,IP AK(7 , 126 » 

• • • 

CLEAR  OUT  J-Th  R^W  IPAK 

• • • 

00  1 K=l,7 
IPAK(K,JI=. 

PAL<  1PR(129)  INTO  IP AKt K, J> ,K=1,7 

• • • 

Ou  2 N=l,12« 
ilN=  (N-  1)  /2  . 

Hi.IsMN*! 

lSH*FT  = 2:-(  't-EC’Hr  I 
KSHIFT=3*1 jHiFT 

IPA<«iNl,J)  = lPAK(Ki,i,  J),oK.  SHIFT  lIPR(N),KaHlFT) 

RETURN 

ENu 


f 


wrs  Uji  i 

con  noos 


PROGRAM  XMARK 


SU6«3UTIN£  iUKFIT 


72/T‘ 


Oi>Tsi 


SUSkOUTInF  CJ-\FII  (FH,F3,FF,ufcL,iNrP£» 


OETcKMINE.  NFW  TmIAl  t//LUt. 

• • • 


A=JFF»FM-2. 3*F0I  / Ic.  u •0£L*t'£,L) 
3=lFF-Ft1»/«2.0*0:L> 

IF(A)1,2,3 


POSITIVi  ?A-?a30lA  ( USF  INFllCTICN  FOIKT. 


10 


15 


31 

32 
C 


• • • 

ZIP=-e/( 2. 0*A) 
IF(ZIP»  31»31,32 
INi,f<t=2IP-r.E 
KeruoN 

1nCSl  = ZIP»''.5 
RETURN 


NtUAIIVi  PA9A30LA  t TEST  SLjFc 


20 


IF  lE.LT.O.I  GC  TC  11 


POSITIWi  SLoPt  t USt  (♦»  FORMULA. 


25 


1NCkE= I-3*5QKT (6*6-4. n*A*FO> )/ (2.0*A)*.P 
RE  TURN 


NEGATIVE  SLOFF  ( USE  (-»  FORMULA) 

• • • 


11 


o 0 


INLwt= (-3-S^RT( 3*6-4. " *A *FC ) ) / ( 2 . 0 * A ) ♦ . 5 

return 


STRlGhT  line  INTEFCEPT 


• • • 

lNCRt=-3/'’0 

RETURN 

FNu 


SYMJOtlL  REFERENCE  MAP  (R=l) 


ENTRY  POINTS 
3 CUR-IT 


VARIAULES 


Sn  t r P E 


Kt  LOLATION 


51 

u 

R-AL 

ol 

3 

PPA. 

0 

OEu 

R - AL 

F.  P. 

0 

FM 

REA^ 

3 

FJ 

PEAL 

F.  F. 

0 

FF 

REA. 

0 

INGRE 

INT  EGER 

F.F  . 

62 

ZIP 

PE  Aw 

EXTERNALS 

SQRT 


T YPt 
■<f  A L 


ARGS 

1 LiRRAPy 
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PROGRAM  XMARK 


SuatJUTlNc.  SIFT 


Tint*  0?T  = 1 


SUBROUTIN'  sift  (NPIXN «NPIX0I 

SEARCH  FOR  TARGET  mS  A 3 X 3 ARRAY  OF  K3T  (3. OR. 4) 
OINENSION  NPIXN(7«12a) . NPI XO ( 7 , 1 28 ) 

LOOP  OVER  tNTIRE  PICTURE,  LXTRACT  I,J-TH  PIXEL, ?USH 

• • • 

00  1 J>1,128 

ZERO  OUT  STACK 
• • • 

Kl=3 

K2=a 

K3=0 

00  1 1*1,128 
K1=K2 
K2=K3 
K3=0 

call  fetch  (NPixa,!, j,kp) 

IF  (KP.LT.3.0R.KP.GT.4)K3*1 
TEST  STACK 

IF  (K1^K2»K3  .EQ.3)  f O TD  7 
CONTINUE 
RETURN 
J2= J*1 
J3= J*2 
11=1-2 
12=1-1 

CALL  FETCH  (NPIX0,I1, J2,KP) 
lF(KP.Ea.3.0R.KP.Ea.4)G0  TO  1 
CALL  FETCH  I NPIXO, 12,  J2, KP) 

IFfKP.E3.S.QR.KP.EQ.4)G0  TO  1 
call  fetch  tFPlXO,!  ,J2,KP) 

IF(KP.E3.3.0K.KP.EQ,4) GO  TO  1 
CALL  FETCH  (NPIXO, II, J3,KP) 

IF(KP.EQ.3.0P.KP.EQ.4IG0  TO  1 
CALL  FETCH  ( NPIXO , 1 2, J 3, KP) 

IF(KP.EQ.3.0R.KP.EQ.4)G0  TO  1 
call  FETCH  (NPIXO, I ,J3,KP) 

IF(KP.E(3.S.UR.KP.LQ,4)G0  TO  1 
• • • 

TARGET  FOUND  - IHSERT  CROSS 

ICR0SS«I-1 

JCR0SS=J»1 

00  88  Ksl,7 

NPIXN(K, JCROSS)=0 

HASK«77777 7777777777777708 

ICOL=(ICR0SS-l)/20 

ICOLl=i:OL*l 

NRES=ICR3SS-20*ICOL 


NSHIFT*(20-NRES)*3 
NASK=SHIFT(NASK, NSHIFT) 
00  99  K* 1,128 
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,^ls  PAGV  rc  U ., 

*•““  cof  I f“aiCA«., 


3JB<3UTINr  SIFT 


PROGRAM  XMARK 

T3/7’-  0PT  = 1 


B9 

C. 


NPaXNJ  luOL  1 ,K)  =NFi  XN(ICOLl«KI.ANC,*ieSK 


b c 


• • • 

PKINI  9£SULI 


C 

i; 


65 


P-vINT  17  ,IC90SS.  JCPoSS 

FuivHAT  ( 23rt3******TARGtT  AT  ( I , J)  = ( iJ,  IH,  I 3 , 7H  J 

rvbT  UPN 

ENC 


V) 


SU9<3UTINl  P'TOm 


7377*.  0PT  = 1 


SOBkOUTINE  fetch  (1PAK,I,J,K) 
JIKENSION  IPAK  <7,126) 


PtACE  I.J-TH  PIXFL  FROM  IPAK  INTO  K 


10 


15 


IF  (I.LT.1.0^.I.GT.l2e )3u  TO  7 
IF  ( J.LT .1 . JR. J.ST.12fl ) SO  TO  7 
Mi<=  (I-l)  /2C 
MNI=MN*1 

NSrtIFT=( 1-2 j*MN) *3 

K=SHiFT(I3AK(HNl,J),NSHIF7).&N0.73 

RETURN 

K=3 

RETURN 

£NL 


SUaxDUTINE  :uRP2  7377*.  0=J=1 


SUoROUTIN£  SOR,2(IFIxi,IF>i.*£,C0lF,bUM) 
OIhe,nsICN  IPIXl(7,l2o),I?IXH7,126) 
iUM=0 . 

00  50  1=1,, 

DO  5i,  J=1,125 
00  5C  K=9,M,3 

N1FIX1=3HIFT(IPIX1(I,J),K-1).AN0.76 
NIFIX2=SMlFT(iOIx2(I,JI,K-l).AN3.73 
CIFF=NI>IX1-NIPIX2 
5.  SUH=SUM»A SSCUIFF) 

C0bF  = SU'171  1 Hfabf-  . 

cot  F=i.:-:u'-F 

RE  TUPtJ 
END 


PROGRAM  DISPLAY 


01SPL4Y  74/74  0^1  = 1 


1 


7 


13 


15 


23 


25 


33 


35 


•*3 


45 


63 


55 


C 

c 

c 


c 

c 

r. 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 


P^OGKft  ^ 01 SPLAY( INPUT, OUTPUT,! #P El .TftPf 3 ,T ftPE 5 1 = CUTPU T 3 
or NtNSION  SHADES f 8), PLiNi (120 ) ,NP1Y( 80 ,643 , TITL ES (7 ) , CHR ( 12 » 
COM'iON/BLKl/7  (200  3 

DATA  SHADrS/l«4,lRP,lRC,lR3,lRC,lP2,lRl,lR  / , 

1 TITLtS/lOHSUBPI Y ,10H0ISPLAY  , 

1 IJH  ,10HFfiAHE  NO., 

1 lOHCHftK ACTtRS, 1 JHUSED  , 

1 10H{  ia3N/P43/ 

CALL  HOOESG(7,33 
CALL  PAG£G ( 7,0 , 0, 13 
call  CPTIO (Z,12HN9C  GH  A1153 
CALL  PAGtG(Z,J,0,2) 

CALL  SETS«G(7,19,4095.  ) 

CALL  ScTSHG{Z,2  J,  3071.  3 
« • « « « 

INPUT  SUBPIX  DATA  INTO  NPIX  ARRAY 

* • * 

RcADi 1,623 ( (NPIX ( I, J 3,1=1, 803 ,J=1, 643 
62  FOl- HAT(20203 

***** 

CAKO  INPUT  PARAMETERS 

13  P'C  - END  OF  FILE  PARA3'ETEF(  HOLERITH  CHARACTER  3 
...  IF  PRC  = ~S”  THEN  PROGRAM  EXITS. 

...  IF  PRC  IS  ANY  other  character  THE  PROGRAM 
CONTINUES. 

23  K1  - SUBPIX  C00RD1NATE(  OCTAL  NUMBER  3 

31  KS  - SUBPIX  COOROINATE(  OCTAL  NUMBER  3 

•*3  K3  - SUBPIX  COORDINATt(  OCTAL  NUMBER  3 

53  NUM  - FRAME  NUMB;R(  INTEGER  3 

63  ILY  - FRAME  ORIENTATION  PAPAMETER(  INTEGEf?  3 

...  IF  ILY  = 0 then  DATA  IS  FOR  TOP  HALF 

...  IF  ILY  = 1 THEN  DATA  IS  FOR  BOTTOM  HALF 

73  ITYPE  - SINGLE  OR  OUADRUPLE  PRINT  PARAMETER(  INTEGER  3 
...  IF  ITYPL  = 0 THEN  SINGLE  PRINT 

...  IF  ITYPL  * 1 then  quadruple  PRINT 

***** 

19  READ  23,  PRC, Ki,K2,K3, NUM, ILY, ITYPE 
2..  F0PMAT{A1,302,I2,2I13 

IF (PRC,tQ.lHS3  GO  TO  99 
CALL  CCNVERT(K1,KK13 
CALL  CONVERT (K2,KK23 
CAuL  CONVt  RT(K3,KK33 

CHECK  TO  SEE  IF  PROCESSING  TOP  OP  BOTTOM  SUBPIXS 


IF(IlY3  30,25 

25  CALL  HEAOE R(TITLES,NUM3 
XX=  514. 

YV=2556, 

CALL  PIX0IS(XX,YV,KK1,KK2,KK33 
X=64, 

Y=2^56. 

IF(ITYPE3  27,29 
27  XHARsZE, 

YB AR=2445. 

29  Kl=K14l 
K2=K2+1 
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PROGRAM  DISPLAY 


PROGRAM  DISPLAY  74/74  0PT=1 


6) 


ti 


73 


75 


83 


85 


9/ 


95 


lu3 


1C5 


<3=<3»1 

CALL  Sf;TSMG<Z,45,,75  ) 

CALL  S'^TSHG<Z,52,2,C) 

DO  1 J=l,40 

00  2 1 = 1,  40 

CALL  TLUd  , J,K1  ,I0X,NPIX» 

2 PLIN^<I)=SHA0ES(IDX+1) 

C «««»« 

00  3 1=41,83 

IX=I-4  C 

CALL  TLU(IX,J,K2,I0X,NPIX) 

3 PLiNL(I)  = SHAOES(IOX+H 


DO  4 1=81,120 

IX=  I-8J 

CALL  TLU(IX,J,K3,IDX,NPIX) 

4 PLINEd  )=SHA0ES  (IDX-f  1) 

C 

CALL  NSTUFF(PLINE,CM^) 

CALL  LEGNDG(Z,X,V,120.CHRI 
IFdTYO£)  6,7 

6 CALL  LEGNDG(Z,X,Y0AR,12J,CHR) 

CALL  LEGNDG(Z,XBAR,Y,123,CHR) 
call  LtGNOG{Z,XBAft, YBAF ,12u ,CHfi) 
YBAR=YBAP-23. 

7 Y=V-23. 

1 CONTINUE 

IFdLYJ  43,45 
4C  XX=614. 

YY=466. 

call  P1X0IS(XX,YY,KK1,KK2,KK3) 
CALL  PAGEG ( Z,a, 0, 1) 

45  GO  TO  19 
C 

33  X=b4. 

Y=1536. 

IFdTVDE)  35,36 

35  XBAR=7E, 

Y 9AR  = 1525. 

36  GO  TO  29 
C 

C CLlAN  op  and  exit 

Q ««*»• 

99  call  CRTID  ( Z,15HN90  END  O''  RUh  ) 
CALL  PAGEG(Z,0,a,2) 

CALL  EXITG(Z) 

STOP 

END 


SYMBOLIC  PLFfcPLNCE  ^AP  |R= 1 I 
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PROGRAM  DISPLAY 


L> 


SUBROUTINt  TUU  74/74  OPT=l 


I SUBROUTINE  TLU ( I, J,K, lOX, HP IX 1 

JIHEKSION  NPIX(aO,64) 

LR=I/21 

NShIFT= <I-(LP*20) 1*3 
5 JROWs  ( (J-H*2)  ♦l^LR 

JPIXsNPIXC  JROH, KJ 
lOXsSHIET ( JPIX.MSMIETJ .AH3,  78 
RE  T URN 
END 


subroutim;  header  74/74  opt=i 


1 SUBROUTIN-  H£ADER(HOP«NUH0£R) 

DIMENSION  H0R(7) 

COMHCN/BLKI/Z (200 ) 

CALL  SETSMG(Z,45,1,5) 

3 call  SETSHG(Z»52« 0.0  ) 

X=  ^400. 
y=2]Jl/. 

DO  till  NMMH=1,3 
CALL  LEGN0G(7,X,Y ,10 tHOR (NMMH) » 
IJ  1111  Y=Y-67. 

Y=15S1. 

call  LEGNOG(Z,X,Y,10,HOR(4> i 
Y=Y-67. 

call  NUMB'’G(Z,X,Y, 2, NUMBER) 
li  Y=1328. 

no  2222  NGDH=5,7 
call  LeGHDG(Z,X,Y,10,HOR(NGOH) ) 
2222  Y=Y-67. 

RETURN 

20  END 


SUBROUTINE  PIXOIS  74/74  OPT=l 


I SUBROUTINE  PIXOIS(X, Y,KA,KB,KC) 

COHHCN/BLK1/7(2JO) 
call  SETSMG(Z,45,1,5) 

CALL  SETSmG(Z,52,0.0) 

3 CALL  NUMBPG(Z«X,Y,2,KA) 

x=x»q20. 

CALL  NUMBRG(Z,X,Y,2,K0) 

X=X492n. 

CALL  NUMBRG(Z,X ,Y ,2«KC) 

Ij  RlTUKN 

end 
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I>ROGRAM  DISPLAY 


SUBROUTINE  NSTUFF  74/74  OPT  = l 


1 SUBROUTINE  NSTUFF (PR LI N, CH ARS ) 

DIMENSION  PKLIN(12C) .CHARS (12) 
COMMON/BL<l/7 (2J0 ) 

H = 0 

5 00  5 J<=1.12 

30  5 JL  = 1 ,10 
HsM*'! 

CALL  PUTCZZ(Z,PRLIN(M) , JL . CH AR S ( JK) ) 
5 CONTINUE 

10  RETURN 

END 


SUBROUTINE  CONVERT  74/74  OPT=l 


1 


13 


SUBROUTINE  CONVERTd.IIt 
C 

C SUBROUTINc  CONVERT  TARES  A TRO  DIGIT  OCTAL  NUMBER  AND  CONVERTS 

C THIS  NUMBERd)  INTO  A DECIMAL  NUMBERdI)  THAT  HAS  A DIFFERENT 

C ABSOLUTE  VALUE  BUT  HAS  THE  SAME  DIGITS  SO  WHEN  THE  OUTPUT  IS 

C PRINTED  THt  NUMBER  APPEARS  IN  CCTAL  FORMAT, 

C 

JONE=I. AN0.7B 

JTENs (SHIFT (1,-3) .AND. 7B) *10 

11= JTENRJONE 

RETURN 

END 


SVM80LIC  REFERENCE  HAP  (R=l) 


ENTRY  POINTS 
3 CONVERT 


VARl AB.ES 

0 I 

SN  TYPE 
INTEGER 

RELOCATION 

F.P, 

0 

II 

INTEGER 

13  JONE 

INTEGER 

14 

JTEN 

INTEGER 

INLINE  FUNCTIONS  TYPE  ARGS 

SHIFT  NO  TYOE  2 INTRIN 


STATISTICS 

PROGRAM  LENGTH  IBB 
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